Compact analytical form for a class of three-loop 
vacuum Feynman diagrams 
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O ■ Abstract 
00 

Q\ ; We present compact, fully analytical expressions for singular parts of a class of three- loop 

diagrams which cannot be factorized into lower-loop integrals. As a result of the calculations 
we obtain the analytical expression for the three-loop effective potential of the massive 0(N) 
Ph! (p A model presented recently by J.-M. Chung and B.K.Chung, Phys. Rev.D56, 6508 (1997). 

jjn : PACS numbers: ll.10.Gh 



Though Quantum Field Theory already has a long history and a number of different ap- 
proaches, Feynman diagrams (FD) are still the main source of its dynamical information. Vac- 
uum (bubble) FD (without external momenta) considered here have several points of applica- 
tions. First of all, it is the evaluation of the effective potentials (for the recent study see Jl j 
and references therein) and renormalization group characteristics (/5-functions and anomalous 
dimensions) of quantum field models and specific operators (for example, anomalous dimen- 
sions of operators in the Wilson expansion). The second important place of the applicabil- 
ity is the calculations of various processes (essentially in Standard Model), where there is a 
possibility to neglect most of the masses and momenta and to calculate only some Taylor co- 
efficients of multipoint FD. These coefficients are bubbles having different (sometimes quite 
big) powers of propagators (i.e. indices of propagators). Using different recurrence relations 
0-[j§] it is possible usually to represent these Taylor coefficients as sets of very simple (usu- 
ally one-loop) bubbles and several so-called master integrals, which cannot be factorized into 
sums of lower-loop integrals. These recurrence relations are some particial cases of the relation 
P, |6| for a general n-point (sub)graph with masses of its lines m 1; m 2 , m n , line momenta 
Pi,P2 =Pi~ Pn,Pn =Pi~ Pin and indeces j u j 2 , -Jn, respectively, 



where Ck = pi + ml are the propagators of n-point (sub)graph. The equation ([!]) is based on 
the rule of integration by part |7]] and has been discovered in the process of the calculation of 
propagator- and vertex-type diagrams (in || for the case n = 3) and of n-point diagrams (in 
for the arbitrary n). 

Because the diagram with the index (i + 1) of the propagator q may be represented as the 
derivative (on the mass m^), Eq.(|I]) leads to the differential equations (in principle, to partial 
differential equations) for the initial diagram (having the index i, respectively). This approach 
which is based on the Eq.(|l|) and allows to construct the (differential) relations between di- 
agrams has been named as Differential Equations Method (DEM). For most interested cases 
(where the number of the masses is limited) these partial differential equations may be rep- 
resented through original differential equation 1 , which is usually simpler to analyse. For the 
bubbles, DEM has been applied in 
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In the case of large number of diagrams (for example, in the calculation of various processes 
of Standard Model) it is convenient ||, || to use Eq.(P without representating the diagrams in 
its r.h.s. as derivatives of initial ones. Then Eq.(|l]) gives the connections between different dia- 
grams decreasing essentially the number of complicated integrals. These integrals (i.e. master 
diagrams) may be calculated [T]J with help of DEM as it was discussed above. 



1. In the present article we consider a class of three- loop bubbles (see Fig. 1) having 
two different masses: the masses of thick lines are one-third of the masses of the bold lines. 
The diagrams are interested because they cannot be factorized into lower-loop integrals (see 
discussions in []T2j ) . Moreover, some of them (J (a), K(a), L(a) and M(a)) give contributions 
(see [U, 0) to the three- loop effective potential of the massive 0(N) <y9 4 -model. The diagrams 
have been studied recently in [jnj using the methods developed in the articles |T^, [njQ. The 
following results for their singular parts have been foundP]: 
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lr The example of the direct application of the partial differential equation may be found in 
2 The singular part of M(a), M (b) and M(c) diagrams has been calculated in pQ] using DEM 



part of M(a) has been found very recently by Broadhurst p5| . 

3 Contrary to Q we use the space D = 4 — 2e and sum the terms 7™ (7 
exponents. 
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is Euler constant) and to 
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where the normaliation factor 



oh 9 — 3e 

m ,m 



(4vr) 6V 7I 2 



^2-) ■ exp (3/2C 2 e 2 ) and /I 2 = (4vr/i 2 )e 7 



and ( n are Euler (^-functions. 

The constants A, B and C have been represented jl2] in the form: 



A = f (1,1), B = f(l,3) and C = /(± ±), 



where 
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The puprose of this short letter is to calculate analytically A, B and C constants and, thus, 
to obtain exact results for the singular parts of the FD presented in Fig. 1. 

Before evaluations we would like to stress that the knowledge of the exact values for the 
singular parts of diagrams is very important, because they determine effective potentials and 
renormalization functions of the quantum field models. Moreover, for various physical pro- 
cesses, the regular parts of diagrams may be evaluated numerically (sometimes with rather 
qood quality) but their singular parts should be known analytically because many types of 
them should be canceled in the end of calculations. 



2. To evaluate A = /(l, 1) we represent the r.h.s. of Eq.@ as the sum of two terms A\ 
and A 2 . 

We introduce the new variable s = (1 — x)/2 and represent A 2 as 



A, 



/o 1 — z Jo 3 + s 2 v 4 

The term A% may be rewritten in the form (y = 1 — t) 

Jo Jo y Jo Ji 1 — t 

Changing the order of integration in the r.h.s. of (^), we have 

f 1 ds , / 1 - s 2 \ 



(4) 
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Thus, the evaluation of 



A = 4 



1 ds 



log (1 - s) + log (1 + s) - log 4 



/o 3 + s 2 

is very simple. Integrating by part, we obtain the final result 

"-->•(!)■ 



(5) 



where (7/2(6*) is Clausen integral JL6 
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Repeating above calculations, we have 
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3. We have obtained results for the singular parts of a FD class in closed analytical form. 
These analytical results are important in the calculations of physical processes of Standard 
Model, where many terms may be involved and the verification and the evaluation of the 
singularities is very important problem. Moreover, the constant A is only one numerical factor 
in the recent calculation Jl| of three- loop correction to the effective potential of O(N) (p A model: 
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After above-mentioned calculation of A 
known fully analytically. 



(i.e. Eq.(^|)), this three-loop correction becomes 
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